Asymptotics of Schwarzschild black hole perturbations 



Anil Zcngiiioglu 
Department of Physics, Center for Fundamental Physics, 
and Center for Scientific Computation and Mathematical Modeling, 
University of Maryland, College Park, MD 207^2, US^ 

We study linear gravitational perturbations of Schwarzschild spacetime by solving numerically 
Regge-Wheeler-Zerilli equations in time domain using hyperboloidal surfaces and a compactifying 
radial coordinate. We stress the importance of including the asymptotic region in the computational 
domain in studies of gravitational radiation. The hyperboloidal approach should be helpful in a wide 
range of applications employing black hole perturbation theory. 
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I. INTRODUCTION 

Linear perturbation theory of black holes is a well de- 
veloped theory [H-Q. The most common approach to 
linear perturbations of a non-rotating single black hole is 
to solve the Regge-Whceler-Zcrilli (RWZ) equations for 
a master function built from a tensor spherical harmonic 
decomposition of metric perturbations d, Q . The RWZ 
formalism plays a prominent role in many studies rang- 
ing from linear stability of black holes Q-lyl to extreme 
mass ratio inspirals |7Hl2|. 

An important problem in numerical time domain cal- 
culations of black hole perturbations is the truncation of 
the infinite spatial domain of Schwarzschild time slices to 
a finite domain by the introduction of a timelikc, artificial 
outer boundary. The presence of such a boundary implies 
an initial boundary value problem. As the boundary is 
not part of the physical solution one tries to construct 
transparency boundary conditions so that the numerical 
solution is as close as possible to the physical solution 
without the boundary. 

The construction of transparent boundary conditions 
is a difficult problem, especially in general relativity [isl - 
[iTj . The RWZ formalism plays a central role in cer- 
tain approaches to this problem, for example in Cauchy- 
pcrturbative matching |18l420| , or the implementation of 
Buchman-Sarbach boundary conditions [21- 23|. 

It is not straightforward, however, to construct good 
boundary conditions even for the RWZ equations. A well 
developed method to deal with this difficulty has been 
presented by Lau f23 - [26j . Lau constructs exact radiation 
outer boundary conditions for RWZ equations based on a 
study of time-domain boundary kernels in Schwarzschild 
spacetime (see also the underlying approach developed 
by Alpert, Greengard and Hagstrom in flat spacetime 
|27|). Lau's method makes it possible to perform accu- 
rate, long-time simulations of RWZ equations in a small 
spatial domain, thereby increasing the efficiency of nu- 
merical calculations significantly. 
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On the other hand, truncating the computational do- 
main for studying gravitational radiation complicates the 
access to asymptotics of black hole perturbations con- 
siderably. In studies of gravitational perturbations the 
main interest lies in the radiation signal as observed at 
future null infinity. In simulations on spatially compact 
domains, elaborate methods are applied on timelike sur- 
faces that represent far away observers to extract a signal 
that is as close to the true signal at null infinity as possi- 
ble [28l-l32|. But even when the background is given and 
gauge independent extraction methods can be conceived, 
it is important to have access to the asymptotic domain 
to capture the nonlinear self interaction of the field away 
from the sources, as argued in Sec. IIVBI 

Including null infinity in the computational domain al- 
lows us to calculate the asymptotics of gravitational fields 
and deems unnecessary the problem of constructing radi- 
ation outer boundary conditions. Only recently, the first 
unambiguous waveform from merging black holes could 
be numerically calculated along null infinity using a mix- 
ture of Cauchy and characteristic methods [s^ HJ] . An 
alternative, efficient approach to include null infinity in 
the computational domain is to compa ctify spacelike sur- 
faces that approach future null infinity |35l439l | . Such sur- 
faces are called hyperboloidal as their asymptotic behav- 
ior resembles the asymptotic behavior of standard hyper- 
boloids in Minkowski spacetime. There has been gro wing 
interest into the study of hyperboloidal foliations |40-|4J| 
and numerical solutions to the hyperboloidal initial value 
problem on given backgrounds |45l452l | . 

The hyperboloidal method has been shown to be very 
efficient in studies of gravitational perturbations |5l|. 
The authors of [5l| solve the Bardeen-Press equations 
based on the Newman- Penrose formalism [s^ [s^]- In 
this paper, we extend the hyperboloidal study of grav- 
itational perturbations to include the RWZ formalism. 
A basic motivation for this extension is that the RWZ 
formalism is a more common tool to study gravitational 
perturbations of a single black hole due to its simplicity. 

The extension of the hyperboloidal method to the RWZ 
formalism allows us to compare tail decay rates of gravi- 
tational perturbations within the Bardeen-Press and the 
RWZ formalisms. The comparison reveals insights into 
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the validity of asymptotic formulae on gravitational radi- 
ation and demonstrates the importance of including the 
asymptotic domain in numerical calculations. 

We extend the study [5l| also with respect to the hy- 
pcrboloidal method such that null infinity can be fixed 
at an arbitrary radial coordinate location, and arbitrary 
coordinates can be used near the black hole [4l|, [s^l ■ We 
find that the ap plication of a transition function sug- 
gested in [H, gives a more accurate spacelike match- 
ing method than the one sug gested in jH These 
improvements are essential for dealing with matter terms 
near the black hole within the hyperboloidal approach. 




II. SCHWARZSCHILD SPACETIME IN 
SCRI-FIXING COORDINATES 

The following section is based on [ilj but evades the 
conformal language used in that reference. We present 
the three steps of the hyperboloidal scri-fixing compact- 
ification in Schwarzschild spacetime: choice of a hyper- 
boloidal time function, introduction of a compactifying 
coordinate, and rescaling of metric functions with a suit- 
able conformal factor. 



A. Choice of foliation 

The Schwarzschild metric in standard coordinates (i, r) 
reads 

g ^ ~Fdf +F-^ dr^ da^, with F := 1 - — . 

r 

Here, dcr^ is the standard metric on the unit sphere. The 
mass of the Schwarzschild black hole is denoted by m. 
Hypersurfaces of constant time coordinate all meet at 
the bifurcation sphere near the black hole and at spatial 
infinity in the asymptotic domain (Fig. [1}. The coordi- 
nates are singular at the event horizon r = 2m. There are 
two common ways to deal with this problem. One either 
introduces the tortoise coordinate and applies outgoing 
radiation boundary conditions close to the black hole, 
or one switches to horizon penetrating coordinates and 
applies an excision technique. 

It is instructive to discuss horizon penetrating coordi- 
nates and to compare them with hyperboloidal coordi- 
nates. For both gauge classes, we change the foliation by 
introducing a new time function 

T = t-h{r). (1) 

All foliations discussed in this paper are of the above 
form. The new foliation respects the stationarity of 
Schwarzschild spacetime. The transformed metric can 
be written as 

g = {~a^+^^(3^)dT^ + 2-f^(3dTdr + j^dr^+r^da^, (2) 



Figure 1: Level surfaces of the standard Schwarzschild time 
coordinate meet at the bifurcation sphere and at spatial in- 
finity. 
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Figure 2: Level surfaces of the ingoing Eddington-Finkelstein 
time coordinate Q foliate the future event horizon, but meet 
at spatial infinity in the asymptotic domain. 

where the lapse a, the shift /3, and the spatial metric 
function 7 read 

with H -.^ ^ called the height function derivative. 

A common system of horizon penetrating coordi- 
nates in Schwarzschild spacetime is based on ingoing 
Eddington-Finkelstein surfaces [stI . [58| . The height func- 
tion derivative for an ingoing Eddington-Finkelstein foli- 
ation reads 



The time surfaces foliate the future event horizon (Fig. [5]) 
and all metric components are regular at r = 2m. 

The construction of hyperboloidal coordinates is anal- 
ogous to the construction of horizon penetrating coordi- 
nates. Comparing Figs. [T] and [2] we see that a height func- 
tion with a suitable singular behavior at the horizon ^ 
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Figure 3: CMC- foliation of the extended Schwarzschild space- 
time with K = 0.4, C = 2 and m = 1. The time surfaces foli- 
ate the future event horizon as well as the future null infinity. 



opens the time surfaces at the bifurcation sphere so they 
foliate the future event horizon. Similarly, a height func- 
tion with a suitable singular asymptotic behavior opens 
the time surfaces at spatial infinity so they foliate future 
null infinity. Expanding the asymptotic hyperboloidal 
condition in a Taylor scries gives |4l| 



2 m 



Taylor 



4m 



O 



(5) 



Here, Ct is a free constant. The explicit form of the 
condition depends naturally on the coordinates in which 
it is expressed. We see that the height function to JS]) 
diverges in the asymptotic domain. 

The first discussion of hyperboloidal foliations in gen- 
eral relativity has been made for surfaces with con- 
stant mean curvature (CMC) fs^H^I. Spherically sym- 
metric CMC surfaces have been constructed explicitly 
in Schwarzschild spacetime by Brill, Cavalho and Isen- 
bcrg [gJ and studied in much detail by Malcc and 
O'Murchadha [iO, IH]. The height function derivative 
for such surfaces reads 



J 



CMC 



where J := 



Kr 



C_ 
'72- 



(6) 



The foliation parameters are the mean extrinsic curva- 
ture K and a constant of integration C. We are not 
interested in varying K or C from one slice to another 
jiol |42| . Instead, we drag the hypersurfaces along the 
timelike Killing vector by keeping the parameters fixed. 

The global behavior of CMC-surfaces depends cru- 
cially on the foliation parameters. For numerical applica- 
tions, we would like to have surfaces that approach future 
null infinity asymptotically (for radiation extraction) and 
penetrate the future horizon near the black hole (for ex- 
cision). These requirements can be fulfilled by choosing 
K > and C > 8m^K/3 [H H. Such a slicing of 
Schwarzschild spacetime is depicted in Fig. |31 



The asymptotic behaviour of the height function 
derivative for CMC surfaces © with K > is given 
by 
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Comparing this with ([5|) we see that the parameter Ct 
is related to the mean extrinsic curvature in lowest order 
by Ct = "ijK. 

The height function derivative for ingoing Eddington- 
Finkelstein coordinates, -ffjEP in Q, vanishes asymptoti- 
cally because ingoing Eddington-Finkelstcin surfaces ap- 
proach spatial infinity. This is also true for maximal sur- 
faces with if = 0, which suggests that the mean extrinsic 
curvature can be regarded, in a sense, as a measure of how 
close the CMC surfaces are to being Cauchy or charac- 
teristic It also controls the strength of the spatial 
redshift effect discussed in IIVCI 

Inserting the height function derivative ([6]) into ([3]) 



gives 



= \J r- + F, /3 = -Ja. 



(7) 



We observe an important difference between Cauchy and 
hyperboloidal foliations. The asymptotic behavior of an 
asymptotically flat spacetime written in a Cauchy foli- 
ation implies a ^ 1 and /? ~ as r — >■ oo. The same 
spacetime in a hyperboloidal foliation satisfies (soj 



0(r) and /3-0(r^) as r 



(8) 



We need to take this singular behavior into account when 
introducing a compactifying radial coordinate in the next 
subsection. 



B. Compactification and rescaling 

To map future null infinity to a finite coordinate loca- 
tion, we choose a compactifying coordinate along future 
hyperboloidal surfaces. We define the compactifying co- 
ordinate p by 



r = — , with 17 



f^(p), 



(9) 



so that the set {fi(/o) = 0} corresponds to infinity in 
terms of the coordinate r. The above transformation is 
written in a generality that allows us to use the standard 
coordinate r in any given domain T) by setting 12 (p) = 1 
for p^V. 

We choose such that asymptotically Vl ^ \/r. We 
introduce the rescaling 51 



11 a, 



(10) 



The rescalcd lapse and shift are not, in general, the lapse 
and shift of the conformally rescalcd Schwarzschild space- 
time. The above rescaling is just a mathematical trick 
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to replace the singular lapse and shift by functions that 
are regular at null infinity. For example, for the CMC 
foliation with the height function (jG]) wc get with (jlOp 



a = VJ^ + FQ^, 13 = -J I 



(11) 



where 



J ■.= nj = 



Kp cn^ 

~3 



The regular metric functions a, /3 and the rescaling fac- 
tor are used in the next section to show the regularity 
of the RWZ equations in hyperboloidal scri fixing coor- 
dinates. 



III. HYPERBOLOIDAL COMPACTIFICATION 
FOR THE RWZ EQUATIONS 

The original RWZ formalism assumes a Schwarzschild 
background in standard Schwarzschild coordinates. Wc 
want to solve the RWZ equations with respect to a new 
time coordinate ([1]), therefore we need a covariant ver- 
sion of the RWZ equations. A gauge invariant formula- 
tion of metric perturbations was given by Moncrief [6^ . 
Moncrief 's formulation also relies on Schwarzschild coor- 
dinates. Gerlach and Sengupta presented gauge invari- 
ant odd-parity RWZ equations in an arbitrary spherically 
symmetric coordinate system (63l |63 | . Their formalism 
has been extended by Sarbach and Tiglio [65|. Martel 
and Poisson completed the covariant, gauge invariant ex- 
tension of the RWZ formalism by presenting covariant 
source terms derived from the stress-energy tensor of the 
matter perturbations [6^. We follow the formalism de- 



veloped by Sarbach and Tiglio [65[. The RWZ equations 
for a master function Z reads 

dtZt = CidrZt + C2drZr + c-^Zt + c^Z^ - a^VZ, (12) 

where an index on Z denotes partial differentiation with 
respect to the corresponding coordinate. The coefficients 
read 

ci = 2/3, 

a2 - 72/32 



C2 



1 



C3 = — (79ta — 7/39,.a -I- a/39,.7 — a9t7 -I- 7a9r/3), 
7a 

7'^a 



(13) 



The potential V depends on whether we solve for odd 
(o) or even (e) parity perturbations. The corresponding 
potentials are given by 



y(°) = ^f/(Z + l)-^ 



^(^) _ A2r2((A + 2)r + 6m) + 36TO2(Ar + 2m) 
(Ar + Qrnj'^r^ ' 

where X = {I — + 2). Here, I refers to the angular 
momentum number from the tensor spherical harmonic 
decomposition of the metric perturbations from which 
the master function Z is built [g^. Note that both po- 
tentials fall off as l/r^ asymptotically. 

In this formalism, metric functions may depend both 
on t and r. We do not need the RWZ formalism in such 
generality. The transformation ([T]) has two properties 
that simplify the coefficients (|13p considerably. First, 
it respects the stationarity of Schwarzschild spacetime. 
This implies that the metric functions are independent 
of T. Second, it leaves the relation a7 = 1 invariant. 
Using these properties we get with respect to the new 
time function t 

drZr = CidrZr + C2drZr + C^Zr + Cd^Zr — a^VZ, 

with coefficients 

ci = 2/3, 



2 / 2 

C2 = a [a 



C3 = a^dr 



,2, , P 
C4 = a^dr o? 



(14) 



The next step is to introduce the compactifying coor- 
dinate p given in (j9l) and the rescaled metric functions 
([TU]). The resulting RWZ equation reads 



ci 



C2 



C3 



drZ-Y — —dpZ^ + -j-^dpZp + ~^^r 



L2 



C4 - C2- 



Zp - d?VZ, (15) 



with 



ci = 2/3, 



.2/„-2 W 



C2 = a a - ^ , 



C3 = c?dp ( A 



a- 

2a / ^,2 



^2 

C4 = a^dp \ — — 



V 



y = — , L = n-pdpn. 



(16) 



The barred coefficients are related to the original coeffi- 
cients via 



ci 



C2 C3 



1 



C4 



2dpn 



C2 
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The RWZ potentials are regular due to their asymptotic 
behavior. They read 



1 



1(1 + 1) 



6mr2 



(17) 



{Xp + 6mQ)'^p'^ 

This ends our discussion of hyperboloidal compactifi- 
cation for the RWZ formalism. The regularity of the sys- 
tem (|15p up to and including null infinity is clear from 
the regularity of its coefficients . 



IV. NUMERICS 
A. The computational framework 

The numerical discretization to solve the hyperboloidal 
initial value problem for gravitational perturbations 
is similar to what has been used in |51[. We use the 
method of lines with 4th order Runge-Kutta time inte- 
gration and 8th order spatial finite differencing. We add 
Kreiss-Oliger type artificial dissipation to the evolution 
equation for to suppress numerical high-frequency 
noise [g^l- The convergence properties of the code are 
qualitatively the same as in [5l|. 

We take a simple Gaussian as the initial perturbation 
and set 



Z{0,p) = 0, Zp(0,p) = 0, Zr{0,p)^e 

Here, pc is the center of the Gaussian pulse and a is its 
width. Experiments with different families of initial data 
with similarly fast fall off or compact support deliver the 
same qualitative results. 

One additional feature in our calculation in relation 
to [sij is the arbitrary coordinate location of future null 
infinity. This is accomplished by the following choice of 
the conformal factor 



(19) 



Here, S denotes the coordinate location of Then 
the inverse of the coordinate compactification © reads 

Sr 



" S + r 
The event horizon is located at 

2mS 

For the CMC foliation determined by (|6]) there is a min- 
imal surface inside the event horizon where the coordi- 
nates break down (4^, 13 • We avoid this coordinate sin- 
gularity in our excision technique by locating the inner 



boundary of our computational domain inside and very 
close to the event horizon. The outer boundary is located 
at p — S. 

The in- and outgoing characteristic speeds for the sys- 
tem (|15p are given by 

c± = ^(-^±«'). 

Inserting the metric functions of a CMC foliation (fTTj) 
with the conformal factor we get 



Along null infinity we have 



c+l 



9 



= 0. 



(20) 



The outgoing characteristic speed at null infinity suggests 
that we need to choose a smaller K if we want to calcu- 
late gravitational perturbations on a larger coordinate 
domain with an explicit time integration scheme sub- 
ject to the Courant condition. Equation (PO)) suggests 
that K and S should satisfy K - 1/S. For the CMC 
calculations we set m = 1,C ^ I, K ~ 0.07,5* = 20. 
Fig. |4] shows the coordinate speeds of characteristics for 
these values of parameters. We see that there are no in- 
coming characteristics into the simulation domain, hence 
no boundary conditions are needed. Outgoing radiation 
leaves the spacetime through the outer boundary. 




P 



Figure 4; The dashed and solid curves represent coordinate 
speeds of in- and outgoing characteristics respectively for pa- 
rameters m = 1,C — 1,K = 0.07,5 = 20. There are no 
incoming characteristics into the computational domain. 



B. Quasinormal modes and tail decay rates 

After an initial transient phase that depends on initial 
data, gravitational perturbations go through a quasinor- 
mal ringing phase followed by polynomial decay. In the 
following, we present results for the simulation of odd 
parity (Rcgge- Wheeler) perturbations with the potential 
P7p . Simulations of even parity (Zerilli) perturbations 
deliver similar results. 
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Fig. [5] shows the quasinormal ringing and subsequent 
polynomial decay in a half-logarithmic plot for an ob- 
server at r = 3m and at null infinity. The ringing lasts 
longer for observers closer to the black hole due to the 
weaker tail signal. It has the form 



Z{t) = a e"'^^^ sin(wiT + (p). 



(21) 



Here, wi and uj2 are the mode frequencies, a is the am- 
plitude and if is the phase. We fit the wave signal 
along r = 3m to the above formula using a simple least 
squares method on the interval r € [60m, 120m]. We find 
wi = 0.373664 and u}2 = 0.088952. These numerical val- 
ues are very close to those obtained by Leaver's continued 
fraction method [gl, H^, which read wi = 0.373672 and 
UJ2 = 0.088963. The error in the measured frequencies is 
dominated by fitting accuracy rather than by numerical 
accuracy, similarly as in [sij . 




500 



Figure 5: Quasinormal mode ringing of odd parity gravita- 
tional perturbations as measured along ,J^'^ (top curve) and 
along an observer at r = 3m (bottom curve). 

From the point of view of asymptotics, the tail part 
is more interesting than quasinormal mode ringing, be- 
cause decay rates at null infinity and at finite surfaces 
are different [t^]- In Fig. [5] we plot, for various observers, 
the local power index defined as 



d\n\^{T,p)\ 
din r 



(22) 



for perturbations with 1 = 2. The function pp(r) becomes 
asymptotically in time the exponent of the polynomial 
decay of the solution. 

Gravitational metric perturbations decay with rate 
— (1+2) along null infinity and — (2^-f 3) at finite distances 
from the black hole, where I is the angular momentum 
number of the perturbation. The dominant I = 2 mode 
is expected to decay as —4 at null infinity and as —7 at 
finite distances. This expectation is confirmed in Fig. [6] 
to a great accuracy. 

There is strong evidence that the asymptotic tail de- 
cay rates for gravitational perturbations in the Newman- 
Penrose formalism at null infinity and at finite distances 
are —6 and —7 respectively [111, [7l|. The decay rate of 
metric perturbations agrees with the decay rate of the 



Weyl scalar ip4 near the black hole whereas the corre- 
sponding rates at null infinity differ. This behavior is ex- 
plained by considering the relation between ip4 and met- 
ric perturbations. An asymptotic formula that one uses 
frequently in radiation extraction algorithms relates ■04 
to second time derivatives of metric petrutbations. This 
relation explains the decay rates —6 for ip^ and —4 for Z 
at null infinity. However, the full Chandrasekhar trans- 
formation relating ■04 to metric perturbations includes 
lower order terms that vanish asymptotically (t^ . [TSj . 
These terms dominate near the black hole and therefore 
the rates for ip4 and Z are the same in that region. Note 
that without having access to null infinity, the agree- 
ment between the finite distance decay rates for '04 and 
Z might be misleading. This observation supports the 
importance of including the asymptotic domain in nu- 
merical calculations. 




2000 4000 6000 8000 10000 

Figure 6: The local power index for vacuum perturbations 
of a Schwarzschild black hole. The locations of the observers 
are from top to bottom: null infinity, 2000m, 300m and 3m. 
We used 2000 grid points and quadruple precision for the 
simulation. 



C. Matching 

Most applications of the RWZ formalism in modern 
literature include the treatment of matter terms, such as 
a point particle. Elaborate algorithms are constructed 
to deal with phenomena like discontinuities and shocks 
related to modelling of matter. Many of these algorithms 
depend on a specific coordinate system typically based 
on Cauchy foliations. In this section we show how to use 
such a coordinate system near the black hole while using 
the hypcrboloidal method in an exterior domain. 

The basic idea is to match a truncated Cauchy surface 
near the black hole to a hyperboloidal surface in the exte- 
rior region [4]| . The matching is performed by a suitable 
choice of the height function derivative that determines 
the time transformation ([I]). We set 



H 



m - 1) 

Hr,, 



P < Pi, 

P > Pc- 



(23) 
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Here, Hi describes a Cauchy foliation, describes a 
hyperboloidal foliation and / is a transition function that 
varies between and 1. We set 



/ 



P < Pi, 

Pi< P< Pe, 
P > Pc- 



In 1411 150| we used 



(24) 



(25) 



With the above choice of /t, the function / is smooth 
ai p = Pe but not at p = pi. An everwhere smooth 
transition function has been suggested in [s^ [H^ 



1 



■ tanh 



S { f -K 


P- 




— tan — 




V2 




- pj 




The free parameter q determines the point pii2 at which 
^T^pi^ ) = 1/2 and s determines the slope of fx at pi/2 

On the resulting matched surface we introduce a radial 
coordinate that coincides with the standard coordinate 
near the black hole and compactifies the exterior region. 
We set 



(i-/) + /(i-§) 



p < Pi, 

Pi < p < Pc 
P > Pc- 



(27) 



As an example, we perform a matched surface evolu- 
tion with an ingoing Eddington-Finkelstein foliation in 
the interior and a CMC foliation in the exterior. We set 
m = l,C = 1.5, K = 0.4, S = 20,1 = 10 and use 2000 ra- 
dial grid points. Our time step is very small, At = 0.001, 
due to the large value of K. The reason why we choose 
such high values for I and K is to demonstrate the spatial 
redshift effect along hyperboloidal surfaces discussed be- 
low. The radial domain is p G [1.8, 20] implying that the 
inner radius is an excision surface and the outer radius 
is at null infinity. The inner and outer transition points 
are pi = 9.08 and po = 12.72. The transition region takes 
20% of the radial grid domain, the rest of the domain is 
shared equally by the ingoing Eddington-Finkelstein and 
the CMC foliations. We use the transition function ((25)) 
with g = 1.3 and s = 1.5. Fig. [7] shows a snapshot from 
this evolution taken at time r = 37m. 

The signal along the Cauchy surface varies strongly 
over the grid, whereas the signal in the hyperboloidal part 
has a low frequency due to the redshift effect (Fig. [7]). 
This effect implies that hyperboloidal surfaces are partic- 
ularly suitable for the accurate calculation of signals with 
strong spatial variation, and also that a small number of 
points is sufficient to resolve the signal which makes high 
order numerical methods very efficient in combination 
with the hyperboloidal approach. The spatial redshift 
effect is stronger for larger values of the mean curvature 
K or higher speeds of outgoing characteristics. 




Figure 7: A snapshot from a matched surface evolution. On 
the left of the dashed lines we use an ingoing Eddington- 
Finkelstein foliation while on the right we use a CMC foli- 
ation. The dashed lines denote the area of transition. This 
plot visually demonstrates the spatial redshift property of hy- 
perboloidal evolutions. 



Decay Rates 
-4.0 , 




300 400 500 600 700 800 900 1000 



Figure 8: The decay rates at null infinity (top curve) and at 
r — 3m (bottom curve). For the dashed curves we use the 
transition function (|25|l . for the solid curves we use (|26|l with 
the same numerical setting otherwise. The smooth transition 
function (I26II leads to a more accurate calculation than (I25II. 



The transition function ([25)1 leads to a more accurate 
evolution than ([25]) . To demonstrate the difference, we 
plot in Fig.|8]the tail decay rates for an / = 2 perturbation 
calculated with the two transition functions. The dashed 
line corresponds to the decay rates as measured from an 
evolution with ([^S]) whereas the solid line is calculated 
from an evolution with (j26p using the same parameters 
otherwise. We see that, for the set of parameters we 
used, the decay near the black hole cannot be captured 
accurately by the evolution with the transition function 
(l25l) whereas the evolution with (1261) is accurate. 



V. DISCUSSION 

We performed a hyperboloidal study of gravita- 
tional perturbations of Schwarzschild spacetime using 
the Regge- Wheeler- Zerilli formalism extending an earlier 
work based on the Bardeen-Press formalism [5l|. The 
motivation for the current work comes from the simplic- 
ity and the wide application of the RWZ formalism in 
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studies of phenomena related to gravitational radiation. 

With the hyperboloidal method we can calculate in a 
simple and efficient way gravitational perturbations of 
Schwarzschild spacetime in time domain using the com- 
mon RWZ formalism without introducing an artificial 
outer boundary into the spacetime. On the physical side, 
we gain some insight into the validity of asymptotic for- 
mulae on gravitational radiation. Studies of decay rates 
of gravitational perturbations at finite distances may lead 
to the erroneous conclusion that the Newman-Penrose 
scalar 'ip4 decays with the same rate as metric perturba- 
tions. While the rates are the same at finite distances, 
they are different at null infinity in accordance with the 
asymptotic formula relating ip4 to second time deriva- 
tives of metric perturbations. This observation empha- 
sizes the importance of including the asymptotic domain 
in numerical calculations of gravitational radiation. 

On the hyperboloidal side, our study employs an im- 
proved coordinatization of Schwarzschild spacetime that 
may be useful in the application of our method to prob- 
lems including matter terms near the black hole. The 
generality in coordinatization is important because many 
algorithms handling matter terms rely on jump condi- 
tions or regularization procedures that have been im- 
plemented in certain coordinate systems. In Sec. IIVCI 
we combined horizon penetrating coordinates in the in- 
terior (U) with a CMC foliation in the exterior ^ but 
other choices are possible. One can, for example, use 
Schwarzschild time slices with tortoise coordinates in the 
interior and slices based on the asymptotic hyperboloidal 
condition ([5]) in the exterior. The implementation of the 
hyperboloidal method in applications of the RWZ formal- 
ism should therefore require only minor modifications. In 
our study of the matching, the smooth transition function 
dm applied in [H, [sl] turned out to be superior to the 
simple transition function ([25)) used in [4l|,[53|- We note, 
however, that the implementation of the matching is an 
additional source of error. Considering that wave extrac- 
tion methods are rather reliable when the background 



has been given, the hyperboloidal method may be less 
efficient than the implementation of accurate boundary 
conditions such as those presented in (23 - [26j . The choice 
of method needs to be made on a case by case basis de- 
pending on the accuracy requirements, time scale of the 
simulation and which features of the radiation signal are 
of interest. 

A natural next step in the hyperboloidal approach for 
the RWZ formalism is to handle source terms near the 
black hole. Among further applications of the method are 
Cauchy-perturbative matching with hyperboloidal sur- 
faces in the exterior (l8l - [20j , or the study of second order 
perturbations of Schwarzschild spacetime at null infinity 
[73, [2EI ■ Perturbative results may not be valid even in 
certain weak field situations |76l - l78l |. It is therefore im- 
portant to check such results by performing comparisons 
with nonlinear evolution. 

Eventually, the main problem one would like to tackle 
with scri fixing is the hyperboloidal initial value prob- 
lem for the Einstein equations [zl, H^l • Recently, Rinne 
presented the first successful numerical implementation 
of hyperboloidal scri fixing for the Einstein equations in 
axial symmetry [8l|. It would be interesting to check 
the validity of linear perturbation theory at null infinity 
by comparing results obtained in this paper with Rinne's 
fully nonlinear calculations. 
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